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-^The  objective  of  this  investigation  was  to  simulate  condi- 
tions of  drifting  snow  previously  observed  in  the  field  in  a labora- 
tory hydraulic  flume  using  a sand- water  analog.  Prototype  conditions 
and  model  results  were  evaluated  to  define  modeling  parameters  that 
would  allow  quantitative  correlation  between  actual  drift  conditions 
and  the  model. 

Although  exact  similitude  is  not  feasible,  this  modeling  approach 
has  proven  useful  in  studies  of  snow  conditions  which  defy  theoretical 
prediction  and  would  be  very  time  consuming  or  impossible  to  evaluate 
in  the  field.  Parameters  examined  in  the  model  include  prototype  di- 
mensions, flow  velocity,  geometric  scale  effects,  structural  porosity, 
and  time. 


Past  Work 
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Various  two-phase  systems  have  been  used  in  the  past  to  investi- 
gate the  mechanics  of  blowing  snow,  including  the  sand-water  analog 
used  here.  Theakston  (15,  l6)  took  a qualitative  approach  and  de- 
veloped a model  which  looked  like  blowing  snow  and  produced  drift 
patterns  resembling  those  observed  in  the  field.  While  largely  ignor- 
ing the  quantitative  aspects,  he  was  able  to  use  the  model  as  a prac- 
tical design  tool  to  study  the  effect  of  various  structural  and 
terrain  modifications  in  the  control  of  snow  accumulation. 

Isyumov  (6),  on  the  other  hand,  developed  theoretical  similitude 
relationships  in  a study  of  roof  snow  loads.  Although  he  did  not 
attempt  to  correlate  the  model  results  to  actual  field  data,  he  was 


-1- 


78  06  09  095 


Distribution  statement  a /)\f1  N /)  /) 

Approved  for  public  release;  U U 

Distribution  Unlimited 


jJ) \ 


WUEBBEN 


t 


able  to  compare  quantitatively  the  effects  of  windspeed  and  wind  ori- 
entation, roof  geometry,  snowfall  rate  and  terrain  roughness  within 
the  model. 


Calkins  (l,  2,  3)  also  used  a similitude  approach  in  developing 
models  of  snow  accumulation  around  various  structures.  In  one  study, 
he  used  the  model  to  examine  possible  snow  control  measures  near  air- 
craft instrument  landing  system  antenna  facilities  in  Alaska.  The 
feasible  alternatives  involved  structural , terrain,  and  site  layout 
modifications  that  could  not  be  easily  addressed  in  the  field.  During 
a study  of  a radar  housing  constructed  in  North  Dakota,  Calkins  de- 
veloped quantitative  relationships  within  the  model  to  correlate  the 
horizontal  drift  dimensions  resulting  from  various  combinations  of 
windspeed  and  the  geometric  scale  of  the  model. 

Experimental  Design 

In  designing  a model  of  drifting  snow,  there  are  a variety  of 
factors  to  be  considered.  On  a qualitative  basis,  the  particles 
should  move  by  a combination  of  saltation  and  turbulent  diffusion  and 
produce  a sand  bed  that  looks  like  a snow  field.  While  an  ideal  model 
would  be  geometrically,  kinematically,  and  dynamically  similar,  exact 
similitude  is  not  practical  with  the  sand-water  analog  and  only  the 
more  important  factors  may  be  modeled. 

A basic  assumption  made  is  that  drifting  occurs  in  dry,  cohesion- 
less snow  since  the  sand  used  is  essentially  cohesionless  in  water. 

For  predicting  long-term  snow  conditions,  Isyumov  and  Davenport  (7) 
suggest  an  integrated  approach  in  which  the  physical  modelings  of 
specific  drifting  events  are  combined  with  statistical  descriptions  of 
meteorological  conditions  in  a mathematical,  model. 


The  importance  of  modeling  the  wind  flow  regime  is  generally 
accepted.  Snyder  (10)  suggested  that  complete  modeling  of  a neutral 
stability  atmospheric  boundary  layer  in  general,  terms  would  result 
from  considering  the  following  criteria: 


1)  Geometric  scaling 

2)  Reynolds  number  scaling  (viscous  effects):  VL/v 

3)  Rossby  number  scaling  (Coriolis  effects):  V/Ltl 
k\  Peclet  number  scaling  (heat  transfer):  Lv/k 

5)  Densimetric  Froude  scaling  (gravity):  V/lgLip^-p^/o^)] 


where  V is  flow  velocity,  L is  a characteristic  length,  V is  kinematic 
viscosity,  is  earth's  angular  velocity,  k is  thermal  conductivity. 
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p is  density,  and  g is  acceleration  due  to  gravity.  The  subscripts  s 
and  v denote  particle  and  fluid  characteristics  respectively. 

The  assumption  of  a neutral  stability  atmosphere  should  be  valid 
in  this  model,  since  drifting  normally  occurs  during  periods  of  high 
winds  and  overcast  skies.  For  the  short  distances  considered  here, 
and  assuming  isothermal  conditions  during  drifting,  the  Coriolis  and 
heat  transfer  effects  can  be  neglected. 


Besides  geometric  scaling  considerations,  this  leaves  Reynolds 
and  Froude  similarity  to  be  dealt  with.  Although  it  is  feasible  to 
satisfy  both  of  these  simultaneously,  it  is  often  not  practical  to  do 
so.  Fortunately,  however,  for  fully  developed  turbulent  or  "rough" 
flow,  the  flow  patterns  become  essentially  independent  of  the  Reynolds 
number.  Most  natural  surfaces  are  aerodynamically  rough  even  at  rela- 
tively low  windspeeas.  A widely  used  criterion  for  determining  fully 
rough  flow  was  proposed  by  Sutton  (ll)  using  the  roughness  Reynold's 
number:  Umz  /v  > 2.5  where  U M = shear  velocity,  v - kinematic  vis- 
cosity, a = terrain  roughness  length,  = e/30  (e  is  equivalent  sand 
roughness1? . 


Based  on  the  above  assumptions,  modeling  of  the  earth's  turbulent 
boundary  layer  might  be  approximated  by  simulation  of  the  terrain 
roughness  (6).  This  may  bo  stated  in  terms  of  the  Jensen  criterion: 


L 

m 

L 

P 


(1) 


where  L is  a characteristic  length,  z is  the  terrain  roughness  length, 
and  m and  p denote  model  and  prototype  conditions  respectively.  Since 
the  terrain  roughness  length  z was  not  measured  in  the  field,  it  will 
be  estimated  based  on  measurements  on  similar  terrain.  Using  data 
cited  by  Isyumov  (6),  the  prototype  terrain  roughness  should  be  in  the 
range  of  0.023  m to  0.09  m.  This  corresponds  to  model  roughness 
lengths  in  the  range  of  0.023  cm  < z _ < 0.18  cm  for  geometric  scale 

ratios  of  1:50  and  1:100.  Because  or^the  appreciable  spread  in  the 
estimate  of  z , as  well  as  the  fact  that  it  is  only  an  estimate,  it  is 
not  possible  ?o  model  the  terrain  roughness  accurately.  Indeed,  there 
is  no  clear  basis  for  using  different  model  roughness  lengths  even 
between  geometric  scales  of  1:50  and  1:100. 

Since  the  roughness  length  measured  for  the  bare  flume  bottom  was 
only  0.001  cm,  a system  of  artificial  roughness  was  developed.  By 
randomly  gluing  pieces  of  crushed  rock  to  the  flume  bottom,  the  rough- 
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ness  length  was  Increased  to  approximately  0.08U  cm.  The  resulting 
boundary  layer  thickness  6 in  the  model,  based  on  velocity  profiles, 
was  estimated  to  be  8 to  10  cm  deep.  This  is  sufficient  to  Just 
envelop  the  tallest  model  (7.6  cm). 

Evaluating  the  roughness  Reynolds  number  for  the  lowest  flow 
velocity  used  in  the  model  (mean  velocity  7^  ■ 28.1  cm/sec)  yields 
U*z  /v  **  18.1,  which  is  within  the  rough  flow  regime  defining  Reynolds 
number  independence. 

Now  one  of  the  major  problems  in  using  the  sand-water  analog  must 
be  considered.  According  to  Mellor  (9),  blown  snow  usually  consists 
of  roughly  equidimensional  grains  about  0.01  cm  in  diameter.  It  is 
apparent  that  it  would  be  difficult  to  maintain  geometric  similarity 
for  any  significant  scale  reduction.  Strict  scaling  would  result  in 
the  very  small  model  particles  being  transported  in  a different  flow 
regime  unless  Reynolds  similarity  is  achieved.  Therefore,  in  order  to 
work  under  the  assumption  of  Reynolds  number  independence,  it  is 
necessary  to  distort  the  model  with  respect  to  particle  size.  In- 
creasing the  particle  size  reduces  the  significance  of  viscous  effects 
relative  to  inertial  effects,  but  provided  that  the  particles  used  are 
not  too  large  and  that  they  still  respond  to  the  more  significant 
aspects  of  the  flow  structure,  it  may  still  be  possible  to. achieve 
adequate  similitude. 

Velocity  scaling  will  now  be  examined.  Since  sand  in  water  is 
being  used  to  simulate  snow  in  air,  density  effects  become  important 
and  the  appropriate  relationship  is  the  densimetric  Froude  number 
cited  in  the  initial  list  of  criteria.  Equating  the  model  and  proto- 
type Froude  numbers  yields  X = 0.0531*  X l/2,  which  results  in  velo- 
city ratios  of  1:130  and  1:1^0  for  geometric  ratios  of  1:50  and  1:100 
respectively.  The  term  X denotes  a ratio  of  model  to  prototype  condi- 
tions. For  the  weighted  seasonal  mean  wind  velocity  of  l6.U  m/sec 
during  drifting  cited  by  Tabler  (12)  field  data,  the  corresponding 
mean  model  velocities  would  be  well  below  the  threshold  of  movement 
for  the  sand  being  used  (7  = 25  cm/sec). 

Since  modeling  based  on  the  the  densimetric  Froude  number  is  not 
possible,  the  velocity  scaling  will  be  developed  based  on  gross  parti- 
cle transport  and  threshold  of  movement  characteristics.  According  to 
Isyumov  (6),  this  can  be  accomplished  by  maintaining  similarity  of 
particle  terminal  fall  velocity  W and  of  critical  shear  velocity 
This  may  be  stated  in  terms  of  model  to  prototype  ratios  as  X^  = X^  = 
X„  . t 
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The  fall  velocity  for  the  sand  used  in  the  model  vas  not  measured, 
but  for  a mean  sediment  diameter  of  0.12  cm  and  a shape  factor  cor- 
responding to  naturally  worn  quartz,  the  terminal  velocity  in  water  at 
20°C  would  be  about  1.2  cm/sec.  Corresponding  information  for  snow 
derived  from  information  given  by  Mellor  (9)  gives  W - 0.5  -*■  1.0  m/sec 
for  a diameter  on  the  order  of  0.1  mm.  This  results  in  a velocity  ratio 
based  on  terminal  fall  velocity  in  the  range  of  1:42  < X <1:83. 

t 

Because  of  the  statistical  nature  of  the  entrainment  process,  it 
is  not  possible  to  pinpoint  a threshold  of  motion,  and  different 
observers  would  likely  estimate  different  thresholds.  But,  based  on 
observation  of  a uniform  bed  of  the  model  sand,  the  mean  flow  velocity 
at  the  initiation  of  particle  motion  is  about  7 ^ = 25  cm/sec,  which 
yields  Um  =1.22  cm/sec.  Taking  the  threshold  of  motion  for  cohesion- 
less snow  to  be  in  the  range  of  21  cm/sec  < U < 52  cm/sec  gives  a 
velocity  ratio  based  on  the  threshold  of  motion  in  the  range  of  1:18  < 

X < 1:43. 
v 


Depending  on  a combination  of  geometric  and  wind  structures  fac- 
tors, the  dominant  form  of  snow  transport  may  be  either  saltation  or 
turbulent  diffusion.  Lacking  prototype  data,  and  because  of  the  uncer- 
tainty contained  in  the  estimates  of  the  parameters  above,  a velocity 
scale  will  be  selected  by  equating  the  threshold  shear  velocity  and 
the  fall  velocity  between  model  and  field  conditions . This  results  in 
the  velocity  ratio  X * 1:42. 

The  Model 


The  objective  of  this  investigation  was  to  simulate  conditions  of 
drifting  snow  and  to  define  modeling  parameters  that  would  allow  quan- 
titative correlation  between  actual  drift  conditions  and  the  model. 

Due  to  the  complexity  of  the  distortion  in  the  model,  it  was  decided 
to  simplify  this  correlation  by  modeling  a snow  fence  which  allows  a 
two-dimensional  analysis. 

The  fence  modeled  is  the  standard  plan  snow  fence  designed  by  the 
Wyoming  Highway  Department  and  installed  along  Interstate  Highway  80 
in  1971.  The  performance  of  this  system  has  been  documented  by  the 
Rocky  Mountain  Forest  and  Range  Experiment  Station  of  the  U.S.  Forest 
Service  (12,  13,  14).  The  fence  is  of  the  horizontal  slat  type  with 
50)t  porosity,  a bottom  gap  of  about  40  cm,  and  a leeward  inclination 
of  15°.  Fence  heights  used  in  the  field  include  1.8,  2.4,  3.2,  and 
3.8  m. 

The  testing  was  conducted  in  a hydraulic  flume  at  CRREL.  The 
overall  length  of  the  flume  is  10.7  m,  with  a working  length  of  7.3  m. 


*V..  = 50  an/sec,  F = 2.8  kg/min. 
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Ideally,  mass  flux  data  would  be  used  to  correlate  the  model 
results.  Originally  it  was  intended  to  make  measurements  of  sediment 
flux  and  concentration  profiles  using  an  electro-optical  probe,  but 
due  to  the  operational  characteristics  of  the  system,  the  results  were 
unsatisfactory.  Because  this  information  was  lacking,  the  model  was 
compared  internally  based  on  an  open  field  (no  fence  in  place)  accumu- 
lation rate  for  various  combinations  of  flow  velocity  and  sand  feed 
rates.  This  accumulation  rate  is  defined  as  the  average  sand  depth  D 
divided  by  elapsed  time. 

The  basis  for  using  this  open  field  accumulation  rate  dD  /dt  is 
the  assumption  that  the  system  approaches  a steady  state.  By  feeding 
in  enough  sand  at  the  headbox  so  that  the  sediment  carrying  capacity 
of  the  flow  is  substantially  exceeded,  net  deposition  occurs  in  the 
flume.  If  a steady  state  of  both  horizontal  flux  and  deposition  rate 
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is  at  least  approached,  deposition  rates  should  remain  more  or  less 
constant  with  respect  to  distance  along  the  flume.  These  accumulation 
rates  were  calibrated  by  running  the  system  at  various  combinations  of 
sand  feed  rate  and  water  velocity.  Keeping  the  inexactness  of  these 
relations  in  mind,  a comparison  of  drift  dimensions  based  on  the 
reference  accumulation  rate  dD  /dt  gives  a relative  indication  of 
fence  efficiency. 

Time  scaling  in  the  model  might  also  be  best  determined  on  the 
basis  of  mass  flux,  but  this  information  was  lacking  in  both  the  model 
and  prototype.  Alternatively  defining  the  deposition  rate  ratio  as 
D^/dt  * *t/V  * xv  yields  prototype  precipitation  rates  in  the 
range  from  33  to  150  cm/hr.  These  model  accumulation  rates  correspond 
to  abnormally  high  values  of  prototype  precipitation  rates  since  most 
snowfall  rates  in  the  U.S.A.  would  probably  be  under  6 cm/hr. 

Moderate  distortion  of  the  precipitation  rate  was  found  to  have 
no  appreciable  effect  on  the  model,  and  even  some  of  the  higher  model 
rates  used  here  appeared  to  cause  only  slightly  more  uniform  depth 
drifts.  As  a result,  it  should  be  possible  to  reduce  the  time  scale 
through  the  use  of  a distorted  precipitation  rate.  If  the  rate  is 
increased  too  greatly,  however,  drift  detail  may  be  obscured  and 
eventually  may  result  in  burying  the  model. 

Another  basis  for  model  comparison  is  the  saturation  value  of 
drift  formation.  The  point  at  which  a fence  becomes  saturated  is  sub- 
jective since  it  is  approached  asymptotically,  and  in  a feed-type 
flume  with  net  deposition  occurring  the  drift  never  really  stops  accu- 
mulating. While  it  appears  that  there  is  a vertical  limit  to  drift 
growth  of  about  1.2  H , sand  continues  to  be  deposited  in  the  tail  por- 
tion of  the  drift.  As  a result,  saturation  in  the  model  was  defined 
as  the  point  when  the  fence  began  to  submerge  rapidly.  Ordinarily  the 
fence  would  remain  relatively  clear  through  most  of  the  run  and  then 
the  bottom  gap  would  plug  and  the  fence  would  rapidly  submerge . 

Whether  this  is  the  best  definition  of  model  saturation  or  not,  it 
unquestionably  signalled  the  end  of  a run. 

Results 


As  a first  step  in  evaluation  of  the  model,  the  results  of  the 
various  model  runs  will  be  correlated  internally,  and  in  addition 
qualitative  comparisons  will  be  made  to  conditions  observed  in  the 
field.  Later,  the  quantitative  correlation  of  the  model  results  to 
prototype  data  will  be  examined. 
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Velocity 

The  data  from  a series  of  runs  using  the  2.U-m  at  1:100  scaxe 
model  demonstrate  the  influence  of  velocity  on  snowdrift  formation. 
Figure  1 gives  a plot  of  dimensionless  area  A/H2  versus  the  reference 
sand  accumulation  depth  D IE.  Since  this  relationship  reflects  a form 
of  fence  efficiency,  vhenrdifferentiated  with  time  it  indicates  that 
the  trap  efficiency  e fence  increases  with  velocity  in  the  range 

of  30  cm/sec  < < 5v  ma/sec,  but  by  = 60  cm/sec  the  trend  re- 

verses. An  examination  of  the  data  for  drift  length  and  maximum  depth 
shows  a similar  trend. 


Figure  1.  Influence  of  velocity  on  fence  efficiency. 


Examining  the  dimensions  of  the  drifts  for  equivalent  areas  A/H2, 
however,  shows  that  as  velocity  increases,  the  drift  increases  in 
length  and  becomes  shallower.  At  * 60  cm/sec,  the  drift  becomes 
both  shorter  and  shallower  by  approaching  a more  uniform  depth.  The 
dimensions  of  the  leeward  drift  at  saturation  for  various  flow  rates 
are  given  in  Table  2.  The  term  denotes  the  distance  from  the 

base  of  the  fence  to  the  point  ofmaximum  drift  depth  and 

denotes  the  distance  from  the  base  of  the  fence  to  the  tail  of  the 
drift. 
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Table  2. 


Variation  in  drift  dimensions  with  flow  velocity 
(saturation  values).* 


M 


w* 


W* 


w* 


A/H 2 


30 

1.19 

2.5 

8.7 

6 

U0 

1.19 

3.0 

2k.  0 

10 

50 

1.19 

3.55 

2U.2 

16 

60 

0.93 

h.O 

26.0 

11 

*H  =2.1 


m. 


1:100. 


Finney  (k,5)  in  his  wind  tunnel  studies,  concluded  that  the 
length  of  the  eddy  area  produced  by  a vertical  barrier  varies  directly 
with  fence  height  but  is  independent  of  wind  velocity.  However,  he 
noted  that  the  distance  from  the  fence  to  the  position  of  maximum 
drift  depth  is  related  to  windspeed.  Kungurtsev  (8),  on  the  other 
hand,  claimed  evidence  of  an  increase  in  a drift  length  with  windspeed, 
and  further  postulated  a linear  relationship  between  drift  length  and 
velocity.  Finally,  Mellor  (9)  cited  data  by  Nokkentved  that  indicated 
a shortening  of  the  drift  as  velocity  increased. 


The  data  from  the  present  study  indicate  that,  for  equal  drift 
volumes,  the  length  of  the  drift  increases  as  velocity  increases  from 
the  threshold  of  motion,  reaches  a somewhat  stable  value  in  the  range 
of  ho  to  50  cm/sec,  and  begins  to  decrease  by  60  cm/sec.  This  would 
correspond  to  the  prototype  drift  length  increasing  with  wind  velocity 


up  to  Vy  * IT  ra/sec  and  starting  to  decrease  with  increasing  velocity 


somewhere  in  the  range  of  21  to  25  m/sec.  Thus,  all  three  trends  men- 
tioned in  the  literature  were  observed  in  the  model  depending  on  the 
velocity  range. 


M 


50  cm/sec,  which  are 


Figure  1 also  shows  that  two  runs  at 
identical  except  for  sand* feed  rate,  agree  very  well.  This  tends  to 
support  the  contention  that  precipitation  rate  is  not  a dominant 
factor  controlling  the  characteristics  of  drift  development,  an  assump- 
tion critical  to  the  use  of  an  abnormally  high  precipitation  rate  ratio 
to  reduce  the  time  ratio.  Data  from  the  runs  using  the  3.8-m  at  * 
1:100  fence  also  support  this  contention.  The  drifts  do  develop  more 
rapidly,  however,  and  there  is  a tendency  towards  more  uniform  depth 
as  the  deposition  rate  increases. 


Fence  Height 


Data  in  which  flow  velocity  and  the  sand  feed  rate  are  held  con- 
stant show  the  expected  trend  of  storage  volume  increasing  with  fence 
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height.  However,  the  3.8-m  at  X^  = 1:50  fence  produced  a larger  drift 
than  the  7.6-m  at  X.  = 1:100  fence  even  though  both  fences  have  the 
same  physical  height.  This  difference  is  attributed  to  the  combined 
effects  of  an  increased  bottom  gap  and  wider  openings  between  slats. 

In  the  field,  porosity  is  often  considered  to  be  the  significant 
variable  rather  than  slat  opening  and  width,  but  because  of  the  parti- 
cle size  distortion  this  may  not  be  true  in  the  model. 

The  maximum  depths  obtained  for  various  fence  heights  became 
stable  before  saturation  and  yielded  the  values  shown  in  Table  3. 

These  data  show  that , a3  fence  height  increases , the  drift  becomes 
relatively  shallower  and  longer,  and  has  its  center  mass  closer  to  the 
fence.  These  effects  were  not  evident  in  the  field  data,  however,  and 
it  is  felt  that  the  problem  lies  in  the  modeling  of  the  boundary  layer. 
Although  the  models  were  contained  within  the  logarithmic  portion  of 
the  velocity  profile,  the  boundary  layer  thickness  was  not  modeled, 
and  the  height  of  the  largest  model  fences  used  approaches  the  thick- 
ness of  the  boundary  layer. 


Table  3.  Variation  in  drift  dimensions  with  fence  height 
(saturation  values).* 


a 

p 

(m) 

a 

m 

(cm) 

XL 

Dmx.'H 

WB 

•YENl/ff 

2.4 

2.4 

1:100 

1.19 

3.55 

24.2 

3.8 

3.8 

1.14 

2.7 

27-8 

7.6 

7.6 

0.77 

1.9 

33.4 

2.4 

4.8 

1:50 

1.07 

3-3 

18.7 

3.8 

7.6 

0.64 

1.6 

30.5 

* 7M  3 50  cm/sec,  Fs  = 2.8  kg/min. 


Porosity 

Another  parameter  to  be  examined  is  fence  porosity.  Although  it 
is  possible  to  geometrically  model  the  fence  porosity,  the  distortion 
present  in  the  model  makes  correlation  impossible  with  the  data  avail- 
able. 


Runs  were  made,  however,  to  determine  the  model's  relative  re- 
sponse to  porosity.  Fences  were  constructed  at  geometric  scales  of 
1:100  representing  a 2.4-m  fence  having  densities  of  50,  75,  and  100£. 
In  addition,  a model  was  constructed  using  square  doweling  to  investi- 
gate the  effect  of  slat  shape. 
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There  is  sl  clear  relationship  between  volumetric  accumulation  and 
.orosSv  ll  leaser  porosities  less  than  50*  as  studied  here  the 

sr 

though  its  geometric  shadow  porosity  is  37*  ana  its  por  a 
inclined  is  50%.  This  indicates  that  that  shape  of 
significant. 

f .‘ss?  ™ 

position  of  maximum  depth  varies  greariy.  n-L’,  observed  in 
duced  the  longest  drift;  this  is  contrary  to  conditions  observed 

the  field. 

Table  U.  Variation  of  drift  dimensions  with  fence  porosity* 

( saturation  values ) . 


Model  Correlation 

the  data  discussed  in  the  previous  section,  it  was  possible 
From  the  da  would  predict  the  unit  storage  volume 

to  develop  a relationship  that  P tiae  based  on  the  fence 

(or  cross-sectional  area  of  the  drifwt  ^ ^ accomplished  by 

height,  flow  velocity  and  sand  feed  r ; portions  of  the  A/E1  vs 

assuming  a linear  ^ff^Sier  d«ermlned  that  the 

DJH  plots  found  in  Figure  1.  It  *»  velocity  was  essentially 

these  two  trends  yields  the  following 

relationship : 

0.162 


D t 

= JL-  x 10 


fl2 


tH 


0.027  + 


- 2.1* 


(2) 


Equation  2 was  evaluated  by  comp. wring  i 
predicted  by  the  equation.  Evaluation  o ^a^a  1Q<  for  ^ost  all 

data^except  of  A/E2.  However,  in  view  of 


II 


WUEBBEN 


the  inherent  inaccuracy  of  using  an  open  field  deposition  rate  as  an 
indication  of  sand  flux,  this  relationship  appears  reasonable. 

It  was  also  possible  to  predict  the  variation  in  the  distance 
from  the  fence  to  the  point  of  maximum  depth  of  drift.  The  relation- 
ship appears  to  be  linear,  at  least  in  the  range  of  30  < 7„  < 60  cm/ 
sec.  Because  of  sparse  data,  the  variation  with  fence  height  was  not 
evaluated,  but  for  the  2.k-m  fence  at  = 1:100,  the  variation  Xyj^ 
with  mean  velocity  may  be  described  as 


*MAX 

■ a 

H 

where  7„  is  in  units  of  cm/sec. 
M 


0.05U  v + 0.36 

M 


(3) 


The  maximum  drift  depth  (£„  „/.¥)  varied  with  fence  height,  but  no 
definable  relationship  was  evident.  The  shortest  fences,  however, 
showed  more  consistent  drift  depths  and,  therefore,  velocity  inde- 
pendence. Although  the  models  were  maintained  within  the  logarithmic 
portion  of  the  velocity  profiles,  the  taller  fences  were  nearly  as 
tall  as  the  boundary  layer  thickness.  The  shorter  fences,  which  per- 
formed better,  more  nearly  approached  the  prototype  boundary  layer 
proportions.  Table  2 listed  the  relative  dimensions  for  the  2.U-m  at 
\T  = 1:100  fence,  which  will  be  used  in  further  analysis  since  it  is 
tfie  best  documented  fence  and  apparently  gave  the  best  results. 


An  examination  of  Table  2 shows  that  the  results  of  the  2.U-m 
fence  are  quite  similar  to  the  prototype  data,  except  that  the  cross- 
sectional  area  of  the  drift  is  slightly  smaller  and  closer  to  the 
fence  than  that  of  the  prototype*  Based  on  the  data  with  the  porosity 
runs  and  considering  the  effect  of  constructing  the  models  from  round 
doweling,  it  appears  that  the  2.U-m  fence  at  a geometric  scale  of 
1:100  with  a nominal  porosity  of  25?  best  represents  the  field  data. 
This  is  dubious  supposition  at  best,  however,  since  there  are  insuffi- 
cient data  to  quantitatively  evaluate  the  influence  of  porosity  in 
this  distorted  model.  Table  5 compares  the  characteristic  prototype 
drift  dimensions  with  those  of  the  model  2.U-m  fences  at  a geometric 
scale  of  1:100  and  porosities  of  25  and  50?. 


Table  5.  Comparison  of  model  and  prototype  data  at  saturation. 


D 

Fence  MAX 

X /H 
MAX' 

xend^h 

A/H2 

Prototype  1.2 

5.7 

2k.  0 

13.0 

Model  25?  porosity  1.2 

5.5 

2k. 0 

17.5 

Model  50?  porosity  1.22 

3.3 

2k.  2 

16.0 
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In  examining  the  formation  of  the  model  drifts  relative  to  the 
field  data,  difficulties  were  encountered,  since  useful  flux  measure- 
ments were  not  obtained  in  the  field  either.  Although  a flux  gage  was 
installed,  it  was  used  more  to  record  the  occurrence  rather  than  the 
magnitude,  of  drifting.  However,  it  is  interesting  to  note  the  simi- 
larity of  drift  formation  in  the  model  and  prototype.  In  both  cases 
the  drift  approaches  its  maximum  depth  rapidly  and  then  the  tail  of 
the  drift  continues  to  build. 

If  the  2.U-m  fence  at  X^  = 1:100  were  taken  as  a reasonable  simu- 
lation, the  corresponding  windspeed  would  be  in  the  range  of  17  m/sec  < 
7^  < 25  m/sec.  This  would  be  a reasonable  velocity  range  since, 
according  to  Tabler's  field  data  (12),  windspeeds  greater  than  17.5 
m/sec  occurred  during  more  than  h5%  of  the  drifting  events,  and  in 
fact  the  mean  windspeed  during  the  drifting  was  l6.U  m/sec.  This 
suggests  the  concept  of  a "significant  wind"  in  modeling,  that  is,  a 
windspeed  that  is  high  enough  to  form  the  snowdrifts  observed  in  the 
field  and  that  occurs  frequently  enough  to  significantly  influence  the 
drift  development. 


Although  problems  with  the  inadequate  modeling  of  the  atmospheric 
boundary  layer  were  encountered,  it  was  possible  to  find  reasonable 
quantitative  relationships  between  the  results  from  different  model 
conditions,  ana  a tentative  correlation  between  these  results  and  the 
prototype  data.  Specific  conclusions  from  the  study  are: 


1.  The  modeling  of  drifting  snow  using  a distorted  sand-water 
analogy  appears  useful  on  a practical  basis.  Because  of  insufficient 
prototype  and  sand  flux  data,  however,  some  assumptions  made  are 
tenuous . 


2.  The  simulation  of  the  atmospheric  boundary  layer  is  of  pri- 
mary importance.  This  may  take  the  form  of  modeling  the  terrain  rough- 
ness and  boundary  layer  thickness.  While  accurate  modeling  of  the 
boundary  layer  thickness  may  not  be  possible,  the  model  structure 

should  at  least  be  small  in  comparison  ( H « 6 ). 

m m 

3.  With  geometric  scaling  based  on  the  terrain  roughness  and 
boundary  layer  thickness,  velocity  scaling  based  on  particle  fall 
velocity  and  threshold  of  motion  characteristics  yields  good  correla- 
tion between  model  and  prototype  data. 

U.  Similarity  of  the  precipitation  rate  is  not  essential,  and  an 
abnormally  high  rate  may  be  used  in  the  model  to  reduce  the  time 


■ 
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scaling.  As  the  precipitation  rate  is  increased,  however,  there  is  a 

tendency  to  obscure  detail  of  the  model. 

5.  The  response  of  the  model  to  porosity  and  dowel  shape  is 

significant,  indicating  a need  for  detailed  model  construction. 
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